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The beam-spin asymmetry, E, for the reaction yd —>■ pn has been measured using the CEBAF 
Large Acceptance Spectrometer (CLAS) at the Thomas Jefferson National Accelerator Facility 
(JLab) for six photon-energy bins between 1.1 and 2.3 GeV, and proton angles in the center-of-mass 
frame, Oc.m., between 25° and 160°. These are the first measurements of beam-spin asymmetries at 
dc.m. = 90° for photon-beam energies above 1.6 GeV, and the first measurements for angles other 
than Oc.m. = 90°. The angular and energy dependence of E is expected to aid in the development of 
QCD-based models to understand the mechanisms of deuteron photodisintegration in the transition 
region between hadronic and partonic degrees of freedom, where both effective field theories and 
perturbative QCD cannot make reliable predictions. 


I. INTRODUCTION 

A. Scaling laws in QCD and experimental data 

The process of deuteron photodisintegration, 

7 + d —f n -I- p, (1) 

is especially important for the investigation of the role of 
quarks and gluons in nuclear interactions. This photonu- 
clear reaction is: i) the simplest {A = 2) and ii) well stud¬ 
ied experimentally. During the last 25 years a number of 
experiments have measured its differential cross section 
over a broad range in energy and angle m- There are 
also some data on the recoil proton polarization laiiQ] 
and the single beam-spin asymmetry mm- 

The most remarkable property of the available cross- 
section data is the energy behavior of this photonuclear 
process. At photon energies > 1 GeV and large pro¬ 
ton scattering angles, it was found that d(T/dt(s, Oc.m.) ~ 
where s and t (and u referred to later in the paper) 
are the usual Mandelstam variables denoting the square 
of the center-of-mass energy and the square of the four- 
momentum transfer to the neutron, while Oc.m. is the 
proton scattering angle in the center-of-mass frame (for 
more details and results on the scaling behavior of the dif¬ 
ferential cross section, see Ref. [I3]). Such a behavior is 
predicted by the constituent counting rules (CCR) based 
on the scaling law for hadron wave functions dHHn]. For 
an arbitrary exclusive two-body reaction at large s and t, 
CCR predict a power-law falloff of the production cross 
section at fixed angles: 

dcr/dt ~/i(6»c.m.)/s”“^, (2) 

where n is the total number of elementary fields in the 
initial and final states, while h{0c.m.) depends on details 
of the dynamics of the process. 
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The quark counting rule was originally obtained based 
on dimensional analysis under the assumptions that the 
only scales in the system are momenta and that compos¬ 
ite hadrons can be replaced by point-like constituents 
with zero angular momentum m [IS- Later, these 
counting rules were confirmed within the framework of 
perturbative QCD (pQCD) up to logarithmic factors by 
showing that exclusive two-body reactions at large s and 
t are dominated by quark and gluon subprocesses at 
short distances m- Within this framework, dimensional 
scaling can be justihed only in the high-energy limit, 
t s TO^, where one can neglect the masses, to, of the 
interacting particles. Therefore, one would not expect 
that the CCR will hold in the few-GeV region. How¬ 
ever, an all-order demonstration of the counting rules for 
hard exclusive processes has been shown to arise from 
the correspondence between a string theory in anti-de 
Sitter space and conformal field theories (AdS/CFT) in 
physical space-time The AdS/CFT correspon¬ 

dence [20] leads to an analytical, semi-classical model for 
strongly-coupled QCD, which has scale invariance and 
dimensional counting at short distances, and color con¬ 
finement at large distances. In this model, dimensional 
scaling occurs not only at very large but also at very 
small momentum transfer, Q, to the parton. In the lat¬ 
ter situation, scaling is due to the constancy of the strong 
coupling with Q when Q is very small. The AdS/CFT 
derivation of the scaling laws is particularly interesting 
since it is a non-perturbative derivation, which suggests 
that dimensional scaling is a feature of both perturbative 
and non-perturbative dynamics. Experimental studies 
of nuclear reactions, such as deuteron photodisintegra¬ 
tion, where the overall momentum transfer is distributed 
among many constituents, so that the momentum trans¬ 
fer per parton is small, are needed to test this model. 

An approximate dimensional scaling has been observed 
in many exclusive reactions at sufficiently high energy 
and large momentum transfer (for reviews see Refs. |21F 
I23j). In addition, the low-energy data on deuteron pho¬ 
todisintegration m (as well as charged-pion photopro¬ 
duction PHI^) also demonstrate scaling behavior. To 
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understand the observed energy behavior it is useful to 
look closely at previous claims of agreement between data 
for differential cross sections and the CCR predictions. 
In fact, the scaled 90° center-of-mass pp elastic scatter¬ 
ing data s^°dcr/dt show substantial oscillations about 
the power-law behavior [26H2^ . Such oscillations are 
also seen in np fixed-angle scattering [3DH32]- The old 
data m as well as the newer data from JLab experi¬ 
ment E94-104 on photoproduction of charged pions at 
0c.m. = 90° [Sll |25] also show hints of oscillation about 
the expected scaling. There are hints of scaling be¬ 
havior in yd —>■ d7r° as well [321 [31] in this case). A 
theoretical interpretation of this oscillatory behavior of 
the scaled cross section was attempted by many authors, 
with the more successful interpretations taking into ac¬ 
count the orbital angular momentum of the partons and 
hadron helicity flip, relating this oscillatory behavior to 
spin-dependent effects. 

The experimental investigations of scaling phenomena 
and related spin-dependent effects resulted in significant 
theoretical advances in understanding the role and range 
of applicability of perturbative QCD at low and interme¬ 
diate energies. These studies make it possible to develop 
a number of nonperturbative QCD-based approaches to 
the hadronic dynamics at long distances. The results 
achieved to date provide a strong motivation for further 
investigation of scaling laws and spin effects in photonu- 
clear reactions through studies of polarization observ¬ 
ables. Measurements of the beam-spin asymmetry, E, 
in deuteron photodisintegration, defined as 

^ 2«<=[E±K*±^3T-'P’4±^6t)-F’2*+^'2-+^5*+^’5-] 

^ — f{e) ’ WJ 

where /(0) = p + |T,_|2] , 

and Fi± = (Ap, A„|T|A.y, Ad) are the helicity amplitudes 
of the reaction, as defined in Ref. [35]: 

Fi± = (±i,±i|r|l,l), F2± = (±i,±i|T|l,0), 

E3± = (±i,±i|T|l,-l),F4± = (±i,T||T|l,l), 

Fsi = (±i,T^|T|l,0), Fe± = (±i, T^|T|1,-1). 


This is done using experimentally determined nucleon 
form factors to describe the gluon exchanges within the 
nucleons. It is hoped that the resulting expressions cor¬ 
rectly include much of the missing soft physics and would 
therefore be valid for momentum transfers lower than the 
ones in the original pQCD expressions. The RNA calcu¬ 
lation is only available at O^.m. = 90° and makes no pre¬ 
dictions for the angular dependence of the cross section. 
It also does not include spin-dependent effects and thus 
cannot make predictions for polarization observables. 


B.2. Hard-Rescattering Mechanism (HRM) 

In the HRM model it is assumed that large-angle hard 
breakup of the deuteron is a two-step process [HI [4T| . In 
the first step, the photon knocks out a quark from one 
nucleon. Then the struck quark undergoes hard rescat¬ 
tering with a quark from the other nucleon, thus sharing 
the high momentum of the incoming photon. Due to the 
hard kernel of the quark-interchange interaction, hard 
rescattering is expressed through the helicity amplitudes 
of high-momentum-transfer nucleon-nucleon [NN) scat¬ 
tering. The number of diagrams accounting for all pos¬ 
sible quark interchanges between the outgoing nucleons 
is very large. However, the HRM allows to effectively 
account for this sum based on the observation that in 
the sum of all possible diagrams, the kernel of the hard 
rescattering can be identified with the quark-interchange 
kernel of the hard elastic NN scattering. The latter al¬ 
lows to substitute the sum of the incalculable part of the 
break-up amplitude with the helicity amplitudes of hard 
elastic NN scattering. 

Given the NN helicity amplitudes, the HRM allows to 
calculate the amplitude of yd —>■ pn scattering without 
any free parameters. One important aspect of the model 
is that while the invariant energy that enters in the NN 
amplitude is the same as the energy of the yd and final 
pn systems. 


These amplitudes can give access to important aspects 
of the underlying physics, such as QCD final-state in¬ 
teractions and quark orbital angular momentum in the 
lightest nuclei. The most popular quark-gluon models 
for deuteron photodisintegration are the reduced nuclear 
amplitudes (RNA) model [36l[37], the hard-rescattering 
mechanism (HRM) [551141] . and the quark-gluon string 
model (QGSM) [3I1I33[. 


sn — Sjd — s — AdJ -|- 2MdE-f, (4) 


the invariant momentum transfer tpf that enters in the 
NN amplitude is less than the one corresponding to the 
yd —pn reaction, t = {k^ — pif)^: 


tN « (fc^ + ^^-Piff 



Ml 

4 ’ 


(5) 


B. Theoretical models 
B.l. Reduced Nuelear Amplitudes Model (RNA) 

The idea of RNA was introduced by Brodsky and Hiller 
in order to extend the region of applicability of pQCD 
down to lower momentum transfers by incorporating 
some of the soft physics not described by pQCD [3S1E7]. 


where toat and Md are the nucleon and deuteron masses, 
E-f is the photon energy in the lab frame, and 
and pif are the four-momenta of the incoming pho¬ 
ton, deuteron, and outgoing proton, respectively (uat « 
{k-y + ^ — P 2 f)^)- This stems from the fact that the 
HRM model corresponds to double scattering in which 
each scattering process carries the half of the total t. 
If we introduce the center-of-mass scattering angle for 
the yd —>■ pn reaction, , and the similar angle for 
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pn —>■ pn scattering, , one obtains from Eq.([^: 


B.3. Quark Gluon String Model (QGSM) 


cos{9 


pn \ 
c.m.) 


1 js-Mj) 
2(s - 4 to^) 

Am% - Mj 
2(s - 4m^)' 


(v^ - \/s - 4to^ cos(02D) 

( 6 ) 


This relation demonstrates that the pn rescattering am¬ 
plitudes enter at smaller angles than 02 m ■ example, 
02.m. = 90° corresponds to 022n. ~ 60° in the elastic pn 
amplitude. 

Using the available experimental values of pn scat¬ 
tering amplitudes, the HRM prediction for the differ¬ 
ential cross section, as well as for the angular distribu¬ 
tion of the hard jd —>■ pn reaction, resulted in a rea¬ 
sonable agreement with the data [a Eg. The very same 
approach also allowed to successfully predict the cross- 
section behavior of the hard break-up of the pp pair in 
the j+^He ^ pp + n^slow) reaction 

Recently, calculations of the beam-spin asymmetry 
within HRM have been updated using the progress made 
in describing pn helicity amplitudes [U |4g. This 
progress is based on the recent observation of the symme¬ 
try structure of valence-quark wave function of the nu¬ 
cleon as well as the new parametrization of the pn ampli¬ 
tude in the large 9c.m. region. The new parametrization 
of helicity amplitudes of elastic pn scattering is based on 
the di-quark model of the valence-quark wave function of 
the nucleon in which the relative phases between scalar 
and vector di-quarks are fixed. Figure shows the s^°- 
scaled pn —>■ pn differential cross section as a function of 
invariant variable s (a) and (b) at 922m. = 69° used to 
describe the energy dependence of the 922m. = 90° helicity 
amplitudes. Using these fits, predictions were made for 



FIG. 1. The scaled differential cross section of elastic pn 
scattering as a function of s and E-/. Data are from Refs. [ig 
148) and the curves correspond to the HRM model fits used to 
calculate predictions for S. 


the beam-spin asymmetry at 922n. = 90° (dotted-dashed 
line in Fig. 13). 

Due to the relation of Fq. ([Hi, angles other than 922m. = 
90° correspond to much smaller values of Iat or ujv and 
as a result, for photon energies between 1.1 - 2.3 GeV, 
tM or um are too soft for the HRM to be valid. 


Another approach to the problem of non-perturbative 
parton dynamics is used in the QGSM proposed by 
Kaidalov [13 Eg ■ Spin variables have been included into 
the QGSM in Refs. EIlEg- This model describes the re¬ 
action through the exchange of three valence quarks with 
an arbitrary number of gluon exchanges. The exchanged 
nucleon is replaced by a nucleon Regge trajectory that 
represents the sum of the exchanged resonances. A non¬ 
linear Regge trajectory provides the best description of 
the data. In a general sense, the QGSM is a microscopic 
(nonperturbative) model of Regge phenomenology for the 
analysis of exclusive and inclusive hadron-hadron and 
photon-hadron reactions at the quark level. Originally, 
the QGSM was formulated for the case of small scatter¬ 
ing angles (j.e., low momentum transfers). Later, Kon¬ 
dratyuk et al. extrapolated the QGSM amplitudes to the 
case of large-angle deuteron photodisintegration [HJ Eg . 
The model fixes all but two of its free parameters from 
other processes, and fixes the remaining two using the 
experimental data on the deuteron photodisintegration 
cross section. It provides predictions for the angular dis¬ 
tribution of the differential cross section and is sensitive 
to spin-dependent effects, making predictions for polar¬ 
ization observables. 


C. Experimental status of deuteron 
photodisintegration 

The extensive studies of the differential cross sec¬ 
tions m have shown that the different theoretical mod¬ 
els describe the available cross-section data on the angu¬ 
lar and energy dependence with about the same degree 
of success. 

Prior to the measurement presented here, there were 
only three sets of polarization data for deuteron photo¬ 
disintegration at energies above 1 GeV. The beam-spin 
asymmetry, E, was measured at Yerevan [m Eg; the 
induced proton polarization, Py, and the polarization 
transfers, Cx' and Cz', were measured at JLab [HEg- On 
the theoretical side, two calculations of the spin observ¬ 
ables are available, within the QGSM [Sg and HRM [Sg 
frameworks. The prediction of the QGSM model for the 
longitudinal polarization transfer Cz' is in good qualita¬ 
tive agreement with the measured data, but the model 
makes no prediction for the transverse polarizations Py 
and Cx' due to their sensitivity to the relative phases of 
the helicity amplitudes [ig. In this respect, calculations 
of Cz' are more stable because they do not depend on 
these phases but only on the moduli squared of the he¬ 
licity amplitudes. The HRM model predictions, as deter¬ 
mined by a parametrization of the pn helicity amplitudes, 
are in qualitative agreement with the available data for 
both Cx' and Cz'- 

For the beam-spin asymmetry, E, there are only the 
Yerevan data in the energy range 0.8 - 1.6 GeV and at 
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0c.m. = 90° [UlITl] (see Fig. 13). Unfortunately, the data 
at « 1.4 - 1.6 GeV have large uncertainties and do 
not allow us to constrain the available models. Never¬ 
theless, the Yerevan data indicate that E(90°) might be 
about 0.5 at these energies. In fact, the QGSM is able 
to accommodate a large beam-spin asymmetry of 0.5 at 
Ey « 1.6 GeV and 0c.m. = 90° [S3], while the HRM is 
not able to do so [SU [SSj . 

The data we present here on the beam-spin asymmetry, 
S, were obtained in an experiment that took place at 
JLab. Our results for E cover photon energies between 

1.1 and 2.3 GeV and nearly complete proton center-of- 
mass angles (between 0c.m. = 20° and 0c.m. = 160°). A 
description of the experimental setup is given below. 


II. EXPERIMENTAL SETUP 

Deuteron photodisintegration was studied using the 
CEBAF Large Acceptance Spectrometer (CLAS) [56] . 
which was housed in Hall B at JLab. GLAS provided ef- 



FIG. 2. A three-dimensional view of CLAS showing the torus 
magnet, the three regions of drift chambers (R1-R3), the 
Cerenkov counters (CC), the time-of-flight detector (TOF), 
and the electromagnetic calorimeters (EC). The CLAS refer¬ 
ence frame, also indicated in the here, was defined with the z 
axis along the beam line and the y axis perpendicular to the 
horizontal. Figure taken from | 56| . 


iment) [53], and electromagnetic calorimeters (EG) to 
identify electrons and neutral particles [60] . 

The geometry of CLAS allowed particle identification 
and momentum determination in a large portion of the 
full solid angle. Charged particles with laboratory polar 
angles between 8° — 140° (this range varies depending 
on the target length and position) were tracked over ap¬ 
proximately 83% of the azimuthal angle with 1-mrad po¬ 
lar and 4-mrad azimuthal angular resolutions. A current 
of —1500 A in the torus magnet produced a magnetic 
field that bent negatively-charged particles away from 
the beamline. The charged-particle tracking system pro¬ 
vided momentum resolution of about 0.5%. Real-photon 
experiments made use of a start counter (ST), which was 
composed of 24 scintillator paddles that surrounded the 
target [61] . The start counter was used in the event trig¬ 
ger and to determine the time at which nuclear reactions 
occurred in the target. 

A linearly polarized real-photon beam was produced 
via coherent bremsstrahlung using a 50-/rm thick dia¬ 
mond radiator, which was positioned on a goniometer. 
The photon beam was then strongly collimated to en¬ 
hance the linear polarization. The characteristics of the 
photon energy spectrum, such as the position of the 
coherent peak, and the degree of photon polarization, 
were controlled by the incident electron energy and the 
orientation of the crystal radiator with respect to the 
beam [53] . Electrons that produced bremsstrahlung pho¬ 
tons were analyzed in the Hall-B tagging spectrometer 
(tagger) [63], which consisted of a dipole magnet and 
scintillator hodoscopes. The tagger allowed the determi¬ 
nation of the incident photon energy by identifying the 
hit position of the scattered electron in the hodoscope 
plane. It provided a tagging range between 20% and 
95% of the incident electron-beam energy. The size of 
the scintillator paddles varied such that an energy reso¬ 
lution of about 0.1% of the incident electron-beam energy 
was achieved. The time of the scattered electron in the 
hodoscope plane was also measured with a resolution of 
better than 150 ps and was used to identify the photon 
that initiated the event detected in CLAS [53] . 

The target used in this experiment was a 40-cm-long, 
conically shaped cell, with a radius of 2 cm at its widest 
point, filled with liquid deuterium. The target cell was 
placed such that its downstream end cap was at the cen¬ 
ter of CLAS. 


ficient detection of particles over a large solid angle. Six 
superconducting coils produced a non-uniform toroidal 
magnetic field and divided CLAS into six identical mag¬ 
netic spectrometers (sectors) as shown in Fig. Each 
sector contained three regions of drift chambers (Region 
1- Rl, Region 2 - R2, and Region 3 - R3) that were 
used to track charged particles and reconstruct their mo¬ 
menta m, scintillator counters (TOF) for particle iden¬ 
tification based on time of flight [55], Cerenkov coun¬ 
ters (CC) to identify electrons (not used in this exper- 


III. EVENT SELECTION AND REACTION 
RECONSTRUCTION 

The data used for this study were obtained during the 
CLAS gl3b data-taking period, which was part of the E- 
06-103 experiment [51] and took place from mid-March 
through June 2007. During this period about 30 billion 
triggers were recorded using a linearly-polarized photon 
beam. The photon-polarization vector was rotated be¬ 
tween two orthogonal directions: parallel and perpen- 
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dicular to the horizontal detector mid-plane, referred to 
as Para and Perp, respectively. Data for six nominal 
coherent-edge positions, 200-MeV apart between 1.3 and 
2.3 GeV, were collected. These data were collected using 
8 different incident electron-beam energies as shown in 
Table H 


(GeV) 

Fe (GeV) 

1.3 

3.302, 3.914, 4.192 

1.5 

4.065, 4.475 

1.7 

4.065, 4.748 

1.9 

5.057 

2.1 

5.057, 5.157 

2.3 

5.157 


TABLE I. Different electron beam energy settings used for 
the six nominal coherent-edge positions during gl3b. 

The trigger during gl3b was relatively loose, a single- 
charged-particle trigger, which led to accumulation of 
data for a number of photoproduction reactions. In this 
study, all events with only one positively-charged track 
were analyzed based on the missing-mass technique. Be¬ 
low we give a detailed description of the procedure fol¬ 
lowed to reconstruct the reaction yd —> pn. 


used as cuts. The proton identification cut is indicated 
by the black curves in Fig. The diagonal bands in 



p (GeV/c) 


FIG. 3. (Color online) Afi as a function of p. The lines show 
the ±3 (t cut from the mean applied to identify protons. 

Fig. are formed by events that were assigned the wrong 
mass, mnorn (Le., non-proton events), as well as acciden¬ 
tal events that were due to particles that did not originate 
in the same physics reaction as the trigger particle. 


A. Proton identification 


B. Photon selection 


Proton identification was done by comparing two in¬ 
dependent estimates of the detected particle’s speed (in 
units of the speed of light, c): one, Pmeas, obtained as 
the ratio of the measured path length from the vertex to 
the TOF and the measured time of flight, and the other 
obtained from the measured momentum and an assump¬ 
tion about the particle’s mass {rrinom)- The difference 
between the two independent estimates was constructed 
as 


AI3 — Pmeas 



(7) 


To identify the protons in our sample, rrinom was set to be 
the nominal mass of the proton. Figurej^shows the event 
distribution of Af3 as a function of the particle’s momen¬ 
tum, p. Proton events are clustered around A/3 = 0. 
The proton-identification procedure was refined by ac¬ 
counting for the dependence of A/3 on the momentum 
resolution of the detector. This was done by dividing the 
distribution shown in Fig. I^into 50-MeV/c-wide momen¬ 
tum bins, and then fitting the A/3 distribution for each 
momentum bin to a Gaussian to determine the mean, 
/iA/ 3 , and the standard deviation, cta^s. The momentum 
dependence of p,A /3 ± 3(ta/ 3 was parametrized and used 
as a proton-identification cut. Fits were not performed 
for the ranges p < 0.7 GeV/c and p > 2.0 GeV/c, due 
to poor statistics, and straight-line extrapolations were 


During the gl3b data-taking period 14 electron hits 
on average were recorded in the tagger for each trigger 
(see Fig. 1^. Thus, for every event, there was a sample 
of ^ 14 photons that could have produced the particle 
detected in CL AS. In order to identify the reaction of 


xIO® 



FIG. 4. Tagger hit multiplicity showing on average 14 photons 
as possible candidates for the true photon that initiated the 
event. 

interest and to calculate kinematic variables, the pho¬ 
ton that initiated the reaction must be selected from this 
sample. This was done by studying the time coincidence 





















7 


between the photon and the proton at the event ver¬ 
tex. The photon arrival time at the event vertex, t-y, 
was calculated using electron timing information in the 
tagger hodoscope, whereas the proton vertex time, 
was calculated using timing information from CLAS. The 
coincidence-time distribution. At = tj — ty, between all 
photon candidates and the identified proton is shown in 
Fig-i The peak centered at At = 0 ns contains photon- 
proton coincident events. The small neighboring peaks 
at 2-ns intervals reflect the bunched nature of the inci¬ 
dent electron beam. The photons in these neighboring 
peaks originate from other beam bunches (not the one 
that initiated the reaction) that came during the trig¬ 
ger window. The photon with a coincidence time within 
At = ±1 ns was selected as the photon that produced the 
proton. Events with two or more photons in this coinci¬ 
dence range were removed from further analysis. Overall, 
photons were unambiguously determined in about 78% of 
all single-proton events. 



effects from these sources, we excluded events in which 
the detected particles fell in a region where the CLAS 
acceptance changes rapidly, by applying fiducial cuts. 
To account for the shape of the drift chambers, these 
cuts were determined by studying the polar angle of re¬ 
constructed tracks as a function of the azimuthal angle. 
Figure shows the angular distribution of protons in the 
six sectors of CLAS along with the applied fiducial cuts 
(black curves). Additional cuts to remove inefficient re- 



FIG. 6. (Color online) Angular distribution of protons as 
detected in the six sectors of CLAS. The black curves indicate 
the fiducial cuts applied to remove events that fell in regions 
where the acceptance changes rapidly. 

gions within each sector did not affect our results since 
those inefficiencies were the same for all Para and Perp 
data, and thus canceled out in the ratio. 


FIC. 5. Coincidence time between all reconstructed photons 
and identified protons. The 2-ns bunch structure of the in¬ 
cident beam is evident. The solid vertical lines indicate the 
±l-ns cuts applied to identify the photon that initiated the 
reaction. 


C. Fiducial cuts 

Charged particles often escaped detection or failed 
track reconstruction in regions near the edges of the 
CLAS drift chambers. Typically, particles that hit the 
support frames or the cryostats of the torus magnet failed 
track reconstruction. In addition, the magnetic field close 
to the torus magnet varied rapidly with position and is 
not modeled very accurately. Therefore, particle tracks 
reconstructed in these regions are characterized by large 
systematic uncertainties, which in turn, propagate to a 
large systematic uncertainty in the reconstructed mo¬ 
mentum. Furthermore, the method used to determine 
the beam-spin asymmetry assumes that the detector ac¬ 
ceptance is constant within each kinematic bin, which 
is not true in the edge regions. To reduce systematic 


D. Event vertex cuts 

In order to reduce background contribution due to 
events not originating in the target, we reconstructed and 
constrained the vertex of each event. The event vertex 
was determined using the distance of closest approach be¬ 
tween the proton track and the beamline position. The 
beamline position was determined for each data run us¬ 
ing multi-charged-track events. Figureshows the event 
distribution over the z-component of the vertex and the 
red lines indicate the cuts we applied to select events that 
originate in the target. The contribution of events with 
X- and y-components of the event vertex outside the tar¬ 
get (greater than 2 cm) was negligible, and for this reason 
no cuts on the x and y vertex components were applied. 


E. Kinematic reconstruction and yield extraction 

The magnitude of the momentum of a charged parti¬ 
cle detected in CLAS was initially reconstructed under 
the assumption that the particle moved with constant 
speed throughout the detector. In order to obtain the 
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FIG. 7. z-component of the event vertex. The solid vertical 
lines indicate the cnts applied to select events that originated 
within the target. 


momentum at the vertex, we corrected this initial esti¬ 
mate for the mean energy loss of the particle as it passed 
through the target, the start counter, and the air gap be¬ 
tween the R1 drift chambers and the start counter [65] . 
as well as for the energy loss in the drift chambers, for 
drift chamber misalignments, and for small imperfections 
in the magnetic field map [66] . Corrections to the inci¬ 
dent photon energies were also applied to account for a 
small gravitational sag in the tagger hodoscope m- 

Deuteron photodisintegration events were then identi¬ 
fied using the missing-mass technique. For each event, 
we calculated the missing mass squared, in the re¬ 
action ■jd —>■ pX using four-momentum conservation. In 
this calculation, the deuteron was considered to be at 
rest, the photon was identified as described in Sec.|^ and 
momentum and energy corrections were applied as dis¬ 
cussed above. Figure [^ shows the missing-mass-squared 
distribution of all events that passed the selection cuts 
described in the previous subsections. Deuteron photo¬ 
disintegration events are clustered in the peak centered 
at the nominal neutron mass squared and were selected 
by the application of a ±3 (t cut on this distribution. 

Figure [^ shows that in the ±3cr missing-mass range 
of interest there was a non-negligible amount of back¬ 
ground, in addition to deuteron photodisintegration 
events. This background contained primarily acciden¬ 
tal events, and events from the reactions yd —>■ pp'K~ 
and yd —> pmr^. The background varied from 5% for 
the low photon-energy bins, to about 40% for the high¬ 
est photon-energy bin. The polarization observable of 
interest could be diluted or altered if the background 
was not removed from the deuteron photodisintegration 
sample. To account for this background, a probabilis¬ 
tic weighting method was implemented. This method 
allowed for a signal-background separation on an event- 
by-event basis in a way that preserved all kinematic cor¬ 
relations [5H1IM] by assigning each event with a signal 


FIG. 8. Missing-mass-squared distribution of the reaction 
yd —>■ pX. The peak corresponds to reactions where the miss¬ 
ing particle is a neutron. The dashed vertical line indicates 
the nominal mass squared of the neutron, and the solid lines 
show the ±3cr missing-mass-squared cuts applied to selected 
deuteron photodisintegration events. 


weight factor, Q, or equivalently, a background weight 
factor, 1 — Q. The Q-factors were then used to weight the 
contribution of each event in the ratio of polarized yields, 
R{4‘) (see Sec.|^. A more traditional approach of fitting 
the missing-mass distribution in each kinematic bin with 
the same predetermined functions was also studied, and 
a comparison of the results was used to estimate the sys¬ 
tematic uncertainty associated with the background sub¬ 
traction method (see Sec. VII. The probabilistic event 


weighting, amongst its many advantages, allows a more 
flexible kinematic binning without having to recalculate 
background contributions, and was thus the method of 
choice for background subtraction in this analysis. 

The signal-background kinematic correlations were 
preserved by dynamically binning data in photon energy 
Ej, and proton angles 0c.m. and (f>. The Q-factor of each 
event was determined by htting the missing-mass distri¬ 
bution of the event’s “closest neighbors” with a predeter¬ 
mined function that described the signal and the back¬ 
ground. The fitting was done independently on Para and 
Perp data, for each nominal coherent-edge position, and 
for each incident electron energy (see Tab.|T]). The size of 
the dynamical bin used, which was defined by the num¬ 
ber of “closest neighbors”, was determined by defining a 
metric in the proton angles 0c.m. and (/>, 


dij — 


cos 9i — cos Oj 

2 

(pi pj 

2 


27T 


1 2 


( 8 ) 


For each event, ej, in the data set, the distances to all 
other events in the data set, d^ (j = l,2...n), were 
computed. Then, a predetermined number of events clos¬ 
est to Bi, Nd, was retained. These events were “closest 
neighbors” of e^. In this analysis, 200 events were re¬ 
tained as closest neighbors. The missing-mass-squared 
distribution of the 200 closest neighbors was constructed 
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and fitted with a Gaussian, describing deuteron photo¬ 
disintegration events, and two exponentials, describing 
background events, 

g{m\,A,^i,a)=Ae "( " ) , ( 9 ) 

6(m|,Ai,A2,Si,B2) = (10) 

These signal and background shapes were chosen since 
a fit using these shapes resulted in the best fit of the 
missing-mass-squared distribution. Other background 
shapes were studied (1®‘ through 4*^-order polynomials) 
without much success in adequately fitting the missing- 
mass-squared distributions. Systematic effects related to 
the choice of signal and background shapes were stud¬ 
ied and are presented in Sec. |VI[ The fit parameters, 
A, /i, cr, Ai, ^ 2 , Bi, and il 2 , determined from the 
fit, were used to determine the signal and background 
functions, g(m\) and b{m\), respectively. The missing- 
mass-squared value of event e^, m^j, was used to calcu¬ 
late the strength of signal {Gi) and background {Bi) for 
that event. 



FIG. 9. (Color online) Missing-mass-squared distribution 
of the reaction yd —>■ pX of Para events in the 1.5- 
1.7 GeV photon-energy bin. The yellow histogram in¬ 
dicates background-subtracted deuteron photodisintegration 
events determined using the probabilistic event-weighting 
method [nn EH]. The red histogram indicates background 
events. 


Gi — 

where N is the number of events in the bin. The statis¬ 
tical uncertainty of the extracted yield in any kinematic 

The Q-factor of event was then calculated to be bin is equal to the sum of the squares of the weights 


Q^ 


G, 

Gi + Bi 


( 11 ) 



(13) 


This procedure was repeated for each event, e^, that 
passed the selection criteria discussed above and yielded 
a unique Q-factor, Qi, for each Ci. Different dynamic 
bin widths were studied with consistent results. A dy¬ 
namic bin of 200 closest neighbors was chosen since it 
corresponds to a kinematic bin width similar to that 
used for the extraction of the beam-spin asymmetry. 
Specifically, for events in the lowest photon energy bin 
= 1.1 — 1.3 GeV), the dynamic bin size of 200 clos¬ 
est neighbors corresponded to about 5° in polar and az¬ 
imuthal angles for events produced at forward angles, and 
to about 10° for events produced at backward angles. On 
the other hand, for events in the highest photon-energy 
bin {Ey = 2.1 — 2.3 GeV), the dynamic bin width of 
200 closest neighbors corresponded to about 10° in po¬ 
lar and azimuthal angles for events produced at forward 
angles, and to about 20° for events produced at back¬ 
ward angles. Figurej^shows the result of the background 
subtraction method for Para events in the 1.5-1.7 GeV 
photon-energy bin. Systematic uncertainties associated 
with the Q-factor determination were also taken into ac¬ 
count and are discussed in Sec. M 

The yield of deuteron photodisintegration events Y for 
each kinematic bin was obtained as 

N 

Y = Y,Q,, ( 12 ) 


In summary, the background-subtraction method em¬ 
ployed here successfully separated signal from back¬ 
ground events while preserving all kinematic correlations. 
Different dynamic bin widths were studied with all of 
them yielding consistent results in the observed quanti¬ 
ties. The uncertainties associated with this method are 
well understood and taken into account. 


IV. PHOTON POLARIZATION 

For the determination of the beam-spin asymmetry, 
the degree of photon polarization had to be known. The 
latter was determined using an analytic bremsstrahlung 
calculation. The photon polarization was found by fit¬ 
ting the enhancement distributions with a theoretical 
calculation m of the coherent spectrum. The enhance¬ 
ment distributions were obtained by dividing the coher¬ 
ent photon-energy spectrum by the photon-energy spec¬ 
trum obtained from an amorphous radiator. The ratio, 
or enhancement distribution, rather than the coherent 
photon-energy spectrum, was constructed in order to re¬ 
move systematic effects, such as counter-to-counter effi¬ 
ciency variations in the tagger, from the determination of 
the photon polarization. Parameters that are character¬ 
istic of the gl3b data-taking period such as electron-beam 
energy, beam collimation, beam divergence and angle, as 
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larization throughout the experiment was of the order of 
75%. Details of the procedure to determine the photon 
polarization can be found in Refs. [siin]- 


V. DETERMINATION OF THE BEAM-SPIN 
ASYMMETRY 


The beam-spin asymmetry, S, is related to the differ¬ 
ential cross section of deuteron photodisintegration as: 


da 

dfl 


da 

dQ 


(1 -I- cos[2r]]). 


(14) 


unpol 


This can be derived directly from the definitions of S 
and the polarized cross section using the helicity 
amplitudes. In Eq. (14) is the degree of linear po¬ 
larization of the beam photon and r] is the azimuthal 
angle between the photon polarization vector and the re¬ 
action plane. In the case of Para events t] = (j) and in 
the case of Perp rj = (j) — 90°, where the angle 4> is the 
proton azimuthal angle measured in the CLAS reference 
frame. The determination of the beam-spin asymmetry 
is simplified by constructing the ratio of polarized yields, 
R{(l)) = [Y(<())ll - Y(((.)-L] / [¥{0 + , where 




0 + 


A<p 






= pll- 


± sin[A<j!)] cos[20]^ A. (15) 


FIG. 10. (Color online) Enhancement distribution (upper 
plot) and calculated degree of photon polarization (lower 
plot) for events with the coherent-edge position at 1.5 GeV. 
The enhancement distribution is htted with the analytic 
bremsstrahlung calculation IZQl and the degree of photon po¬ 
larization is calculated (dashed red curve). The degree of 
photon polarization is then corrected for differences between 
the fit and the enhancement, as well as statistical fluctuations, 
using information from data with adjacent coherent-edge po¬ 
sitions (blue curve). More information on the procedure is 
found in Refs. |6^rfT| . 


well as fluctuations of the coherent radiator position and 
angle, the beam-spot size, and multiple scattering, were 
taken into account in the calculation. Figure shows 
the enhancement distribution for the coherent-edge po¬ 
sition at 1.5 GeV fitted with an analytic bremsstrahlung 
calculation (upper plot). The lower plot shows the calcu¬ 
lated photon polarization based on the enhancement fit 
(dashed red curve) and the corrected polarization that 
takes into account residual differences between the fit and 
the enhancement distribution (blue curve). 

The method and the procedure we used to determine 
the degree of photon polarization produce reliable results 
with small systematic uncertainties in regions where the 
enhancement is large. Therefore, only events with pho¬ 
ton energies above Ece — 200 MeV (where Ece is the 
coherent-edge position for the current event), were kept 
for further analysis. The average degree of photon po- 


Equation (15) results from Eq. (14), with the parameter 
FIIA being the incident photon flux, A(j) the 0-bin width 
used to bin the data, and A the detector acceptance that 
is assumed to be constant within the (()-bin (the effect 
of this assumption is investigated in Sec. VI). The no¬ 
tations II and T indicate the orientation of the photon 
polarization vector. 


Substituting the expression in Eq. (15) in the definition 
of the ratio becomes 


R(0) = 


y(.^)ii - y(^)^ 
Y{(f>)\\+Yi(py 
Fr-1 + 


FrPr+Io PV sin[A(j)] 
ffl + l 


cos[2((). - (^o)] 


Fr + 1+ ^g,^«7^ 2PSg 


where Ed = 


F-L ! 


Pr = 


PF ’ 


and P = ° ° . 


, (16) 


The 


parameter cj^o in Eq. (16) accounts for any systematic 
offset of the photon polarization vector from its nominal 
orientation. 

Acceptance effects and any acceptance related system¬ 
atic uncertainties, cancel out in From the fit of the 

polarized-yield ratio, R{(j))^ to the function 


F{^) = 


^ - 1 + cos[2{(j) - D)] 

A + i+"4|^2 Ccos]2(™)d)]’ 


(17) 


the following are determined: 

• A: ratio of Para and Perp fluxes, Ep, 
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• B: ratio of Para and Perp polarizations, Pr, 

• C\ product of average polarization and asymmetry, 

A<p ’ 

• D: offset of the photon polarization vector, (j)Q. 


The fitting was optimized by fixing three of the four pa¬ 
rameters using independent methods. The optimization 
was extensively studied along with any associated sys¬ 
tematic uncertainties m- Specifically, for each photon- 
energy bin, the parameter A was determined from a fit 
to the azimuthal distribution of R{(j>) integrated over all 
proton angles in the reaction of interest. Since the inci¬ 
dent photon flux was constant for a given photon-energy 
bin, a ht to the integrated kinematic bin ensures adequate 
statistics to precisely estimate the photon flux ratio Pr. 
The parameter B was calculated using the degree of pho¬ 
ton polarization on an event-by-event basis (as obtained 
from the procedure described in Sec. IV ). Finally, the pa¬ 
rameter D was obtained from fits to the high-statistics 
single-pion reaction, 7 c? —>■ PsPT^~ - m- The parameter C 



we used Monte-Carlo data where we could control the 
true value of E. Para and Perp (j) distributions were gen¬ 
erated according to Eq. (141 and binned in the exact way 
CLAS data were binned, removing data that fell outside 
the CLAS fiducial regions. The function of Eq. 0 was 
fitted to the generated azimuthal distributions, and E^et 
was obtained from the fit parameter C as Tidet = CjP. 
Then, the correction factor for the variable (()-bin widths 
was determined by fitting Tidet vs Ege„ with a first-order 
polynomial (see Fig. 12). The slope of the htted line gives 
the value of C.F. For the variable ^-bin widths chosen 
for this analysis, the correction factor determined from 
this study is C.F. = 1.0094 ± 5 x 10“®. Further details 



FIG. 12. Beam-spin asymmetries determined from fitting gen¬ 
erated data as a function of the true value Egen for 36 variable 
4> bins (equal to an average bin width /S.(f> = 13.50°). The 
correction factor due to a variable 0-bin width, C.F., is deter¬ 
mined by fitting Edet vs Egen to a first-order polynomial (see 
text for details). The statistical uncertainties of the points 
are smaller than the symbol size. 


FIG. 11. Ratio for a specific kinematic bin (E~f = 1.3 — 
1.5 GeV and = 120.0° — 132.5°) fitted to the function 
F{(j)) (see text) to determine the beam-spin asymmetry for 
that bin. 


can be found in Ref. na¬ 
vi. SYSTEMATIC UNCERTAINTIES 


was then determined by fitting the ratio i?( 0 ) for a bin 
in 9c m and E.y, hxing all other parameters, as shown in 
Fig.0 

It is common in CLAS data analyses that the 
polarized-yield ratio is distributed in (j) bins of variable 
width (smaller width in regions in the middle of each 
CLAS sector and larger width closer to the edges). This 
complicates the determination of the beam-spin asym¬ 
metry since the correction factor^ sin^At/i] ’ takes only a 
single value of A(j). One would intuitively expect that 
in the case of a ht to a variable- 0 -bin-width distribution, 
the correction factor would be some average over all 0 -bin 
widths; in fact, a good approximation of the correction 
factor can be calculated in this way (i.e., C.F. « ). 

To precisely quantify the value of the correction factor. 


Systematic uncertainties due to various sources were 
studied and quantihed. The uncertainties fall into sev¬ 
eral main categories: due to the (/)-bin method used to 
extract E, due to the choice of cuts for event selection, 
due to the choice of background shape, and due to the 
statistical uncertainties of quantities used to calculate E. 
The parameters of these sources were varied within rea¬ 
sonable bounds and the sensitivity of the hnal result was 
checked against this variation. In several studies, gen¬ 
erated data were used to determine the systematic un¬ 
certainties, and in others, the actual experimental data 
were used. A summary of the systematic uncertainties 
is given in Table |ll] indicating whether the source results 
in an absolute or relative uncertainty. This section sum¬ 
marizes the studies performed to estimate the systematic 
uncertainties. 
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Source 

Estimate 

Type 

Variable cf>-hin width 

10"® 

1 

Detector acceptance 

<1% 

2 

Uncertainty of po offset 

10"® 

1 

Uncertainty of Pr 

1% 

2 

Uncertainty of P 

5% 

2 

Uncertainty of Fr 

~0.002 

1 

Fiducial cuts 

~0.014 

1 

Particle ID cuts 

~0.01 

1 

Missing-mass cuts 

~0.01 

1 

Background subtraction method 

~0.012 

1 

Background shape 

~0.01 

1 

Q-factor 

~ 0.02 

1 


TABLE II. Systematic uncertainties contributing to the to¬ 
tal systematic uncertainty of the beam-spin asymmetry. The 
estimated values are averaged over all kinematic bins. The 
type of the source indicates whether the uncertainty is abso¬ 
lute (Type 1) or relative (Type 2). 


A. Variable (ji-hin width 

The uncertainty of the correction factor that accounts 
for the variable (()-bin width has been propagated to the 
uncertainty of the beam-spin asymmetry using generated 
distributions. The value of this uncertainty was found to 
be of the order of 10“® and is negligible compared to all 
other uncertainties m- 


B. Detector acceptance 

Acceptance effects on the determined beam-spin asym¬ 
metry are twofold. Firstly, the 0-bin method used to de¬ 
termine the beam-spin asymmetry relies on the assump¬ 
tion that the acceptance is constant within each (()-bin, 
and thus can be taken out of the integral (see Eq. (Hll)- 
Secondly, the detector acceptances of the Para and Perp 
data are assumed to be identical, and therefore cancel 
out in the ratio R{4>) (Eq. @). 

The effect of a non-constant acceptance within each (j)- 
bin has been investigated using simulated distributions. 
Details of the study can be found in Ref. [72] ■ The sys¬ 
tematic uncertainty of E due to the assumption of con¬ 
stant acceptance within each (/)-bin was found to be less 
than 1%. 

Data from the gl3b data-taking period was collected in 
a way to minimize any difference between the acceptance 
for the Para and Perp settings. This was accomplished 
by changing the photon polarization between Para and 
Perp about every two hours. Studies comparing the ra¬ 
tio of proton yields in adjacent TOF counters between 
Para and Perp data averaged over all runs within each 


coherent-edge setting show that variations in the detector 
acceptance were within the statistical uncertainties. In 
addition, generated data using different acceptances for 
Para and Perp distributions of the size of the experimen¬ 
tal variations showed negligible effects on the determined 
beam-spin asymmetry. The overall uncertainty of E due 
to acceptance effects is thus less than 1%. 

C. fo offset 

The systematic effect due to the uncertainty of the di¬ 
rection of the photon-polarization vector {(pQ offset) was 
investigated using generated data The uncertainty 
of the beam-spin asymmetry, which stems from the un¬ 
certainty of the (j>o offset, was found to be of the order of 
10“® and is negligible compared to all other uncertain¬ 
ties. 


D. Photon polarization 

The effect of the uncertainty of the photon polarization 
on the estimated value of E is twofold. On one hand, the 
uncertainty of the photon polarization propagates into 
the uncertainty of the polarization ratio Pr, which is 
used as a fixed parameter in the fit (Eq. ([izl))- This af¬ 
fects the fit and the uncertainty of the free-fit parameter 
C that is used to determine E. On the other hand, the 
uncertainty of the photon polarization propagates into 
the uncertainty of the average photon polarization P, 
which is used to calculate E from C: E = ^ . An 

independent study determined that the systematic un¬ 
certainty of the photon polarization was 7% |74j. Using 
the value of 7%, the uncertainty of Pr is calculated to be 


APR = PRy'(^^) 

APr = 0.1 X Pr. (18) 

The uncertainty of P is 

AP = ^(A/j|)^ -I- (APj_)^ and thus, 

AP-0.05xP, (19) 

where we used P^ « pH « P. The uncertainty of E due 
to the uncertainty of P is then 5%. 

To estimate the uncertainty of E due to the uncer¬ 
tainty of Pr, generated distributions were produced and 
analyzed m- The study yielded an uncertainty of E of 
less than 1%. 

In general, the uncertainties due to Pr and P are 
highly correlated and should be treated together. How¬ 
ever, it is evident that the uncertainty of the average 
polarization of ~ 5% has a much bigger effect on E than 
the uncertainty in the polarization ratio Pr (< 1%), and 
the former is quoted. 
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E. Incident photon flux 

The systematic uncertainty due to the uncertainty of 
Fji was determined using studies similar to those for the 
(j)Q offset m From these the systematic uncertainty of 
the beam-spin asymmetry was found to be 

= 0.073(20) 
which on average corresponds to an uncertainty of 0.002. 

F. Reaction selection cuts 

The systematic uncertainty due to the choice of fidu¬ 
cial cuts was determined by varying the cuts from their 
nominal values to tighter values (the fiducial ranges of 
9 and (j) were reduced by ^ 3°). The variation of the 
beam-spin asymmetry was found to be on average 0.014 
and we report this value as the systematic uncertainty 
for this source. 

Proton ID cuts were varied between 2a and 3cr using 
experimental data and the determined beam-spin asym¬ 
metries were compared. The variation of E on average 
was 0.01 and we quote this value as the systematic un¬ 
certainty. 

The missing-mass cuts, which select deuteron photo¬ 
disintegration events, were varied between 2cr and 3cr, 
and the determined beam-spin asymmetries were com¬ 
pared. This study accounts for possible leakage of events 
from background channels due to the non-Gaussian shape 
of the signal. On average the uncertainty of S was found 
to be 0.01. 


G. Background subtraction 


used a dynamic bin width, to results from a binned 
method in which the background subtraction was deter¬ 
mined on a bin-by-bin basis. The study yielded a sys¬ 
tematic uncertainty of 0.012. 

The Q-factor weights have an uncertainty that depend 
on the dynamic bin width as well as on the goodness of 
the ht. Specifically, the calculation of the uncertainty of 
Qi was done by propagating the uncertainties of the fit 
parameters, as 

\" dQi . dQi 

( 21 ) 

jk 

where j and k run over the number of the fit parame¬ 
ters, Pi are the fit parameters, and Cov{i,j) is the co- 
variance matrix determined from the fit. The Q-factor 
uncertainties are highly correlated between events in the 
same kinematic bins due to the method of nearest neigh¬ 
bors. Therefore, the Q-ia.ctor uncertainties for each kine¬ 
matic bin add up, 

^ ' ^Qi ^ ( 22 ) 

i 

where the sum is over the number of events in each 
Oc.m., and (p bin. Through several studies using 
the deuteron photodisintegration events, we were able to 
determine the systematic uncertainty of the beam-spin 
asymmetry, which is due to the uncertainty of the Q- 
factor value, as a function of the statistical uncertainty 
of E, agidi^ 

erg = 0.0093 4-0.176 • cr.tat. (23) 

On average, this corresponds to an uncertainty of 0.02. 
More details on the determination of this uncertainty can 
be found in Ref. [72] . 


The systematic uncertainty associated with the back¬ 
ground subtraction is threefold. Specifically, there is 
a systematic effect associated with the predetermined 
shape of signal and background, which is used in the 
fits and the determination of the Q-factor. In addition, 
there is an uncertainty associated with the background- 
subtraction method. Finally, there is an uncertainty as¬ 
sociated with the Q-factors. This latter uncertainty can 
be determined by propagating the uncertainties of the fit 
parameters to Qi. 

The uncertainty due to the assumption that the signal 
is Gaussian is accounted for in the systematic uncertainty 
associated with the missing-mass cut. The uncertainty 
due to the choice of the background shape was studied 
by comparing results obtained with a linear background 
and with the nominal background (two exponentials - see 
Eq.(lOl). This uncertainty was found to be of the order 
of O.OI. 

The uncertainty associated with the background sub¬ 
traction method itself was studied by comparing results 
from the probabilistic event-weighting method, which 


VII. RESULTS AND DISCUSSION 

Using the data from gl3b, the beam-spin asymme¬ 
try was determined for incident photon energies from 
=1.1 to E.y =2.3 GeV and proton angles between 
9c.m. = 25° and 0c.m. = 160°. In the following sections we 
report the energy and angular dependence of the beam- 
spin asymmetry and associated uncertainties, which were 
determined as outlined in the previous sections. 


A. Energy distributions 

Figure [T^ shows our 0c.m. = 90° data for E compared 
to the available data from Yerevan and to the model pre¬ 
dictions from QGSM and HRM. The precise CLAS data 
are in good agreement with the published Yerevan data 
and increase the kinematic coverage up to photon ener¬ 
gies of 2.3 GeV. The linear energy dependence predicted 
by the QGSM is not confirmed by the data: the model 
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predicts larger asymmetries than the data at lower pho¬ 
ton energies and lower values than the data at higher pho¬ 
ton energies. The updated HRM reproduces the general 
shape of the energy dependence. Especially, it describes 
the increase to higher asymmetries observed in the data 
between photon energies 1.6 and 2.0 GeV. In the HRM, 
this increase stems from features of the pn scattering am¬ 
plitude. However, the model underpredicts the values of 
the asymmetries over the entire energy range. 

Figure [I^ shows the energy dependence of E for four 
different proton center-of-mass angles. The width of the 
photon-energy bin is kept constant at 200 MeV, whereas 
the width of the angular bins varies in an attempt to have 
similar statistical uncertainties. The results indicate pos¬ 
itive asymmetries for angles larger than 50° and negative 
asymmetries for forward-going protons. The 9c.m. = 90° 
result displays the largest asymmetry and suggests a local 
maximum at = 2.0 GeV. 


B. Angular distributions 

Figure [T5| shows the angular dependence of S for the 
six photon-energy bins (200-MeV wide) between =1.1 
and E^ =2.3 GeV, as well as the QGSM and HRM pre¬ 
dictions. The width of the angular bin varies in an at¬ 
tempt to have constant statistical uncertainties. The re¬ 
sults indicate that the beam-spin asymmetry has a local 
minimum at 9c.m. = 90° for the lowest photon-energy 
bins. This minimum evolves to a maximum for the higher 
photon-energy bins. The data exhibit, especially at the 
three lower energies, complex structures. At small an¬ 
gles the observable is negative, but increases with the 
production angle and reaches a positive maximum below 
9c.m. = 90°. Then, it decreases to a positive minimum 
shortly above 9c.m. = 90°, and reaches a second maxi¬ 
mum at large angles. As the photon energy increases, the 
position of the first maximum shifts towards 0c.m. = 90°, 
while the magnitude of the second maximum, observed at 
large angles, continuously decreases and E becomes neg¬ 
ative at these large angles for E^ above 1.7 GeV. While 
the QGSM model predicts similarly complex angular dis¬ 
tributions, there are significant differences between the 
data and the model. The latter predicts positive E at all 
energies and angles, while data show that in some bins 
E is negative. The positive maximum at small angles 
in the model is not confirmed by the data. The positive 
maximum at large angles is confirmed only by the lowest- 
energy data, but at a different angle. The QGSM seems 
to predict well the maxima at 9c.m. = 90°, observed in 
the higher photon energy data. The discrepancies be¬ 
tween the data and the QGSM may be due to resonance 
contributions to the reaction dynamics, which are not 
well described in the model. The HRM model is expected 
to be valid only for 9c.m. = 90° angles, where both tjv 
and un are large. Nevertheless, predictions of the HRM 
model are shown to demonstrate the kinematics where 
the data exclude applicability of the model. 


C. Conclusions 

The work presented here increases significantly the 
kinematical coverage and precision of the available data 
for the beam-spin asymmetry, E, of the reaction yd —> 
pn. State-of-the-art models have limited success in repro¬ 
ducing the details of the data. The fact that the models 
capture only some of the most general features of the 
data suggests that improvements are needed in the de¬ 
tails of the reaction dynamics. For example, a better 
phenomenological input to the HRM for the spin depen¬ 
dence of the elementary pn amplitudes could bring the 
calculation closer to the data. Alternatively, the compar¬ 
ison at 9c.m. = 90° may suggest that the energy range of 
our data are below the full applicability of the model. 
The CLAS data provide stringent constraints that can 
be used in the development of the existing models or 
even aid in the development of new phenomenological 
approaches that attempt to describe the underlying dy¬ 
namics in the transition region from hadronic to partonic 
degrees of freedom. 
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FIG. 13. (Color online) Beam-spin asymmetry as a function of the incident photon energy for 9c.m. = 90°. The red (triangle) 
and cyan (circle) points show the results from Yerevan |11[ I12 | , whereas the blue square points are the results from the present 
work. The solid and dash-dotted lines are the QGSM [53] and HRM | 541I55 | predictions, respectively. The blue band indicates 
the systematic uncertainties of the present CLAS measurements. 
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FIG. 14. (Color online) Beam-spin asymmetry as a function of photon energy for four different proton angles in the center- 
of-mass frame. The angular bin size varies with photon energy. The blue bands indicate the systematic uncertainties of the 
present CL AS measurements. 
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FIG. 15. (Color online) Beam-spin asymmetry as a function of the proton center-of-mass angle, Qc.m., for the six photon-energy 
bins between 1.1 and 2.3 GeV. The solid and dash-dotted lines are the QGSM |53) and HRM [541155) predictions, respectively. 
The blue bands indicate the systematic uncertainties of the present CLAS measurements. 
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